Surface Driven Bulk Reconstruction of Gold Nanorods 
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Molecular dynamic simulations are used to study the heating of a gold nanorod of 2624 atoms. 
We show that roughening of surface {110} facets leads to a shape transformation and structural 
rearrangement of surface and bulk atoms in the rod, in order to reach a more stable configuration. 
Our results show that the stability of the nanorod is governed by the free energetics of the surface 
facets. 
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The properties of metallic clusters of various shapes, on 
the nanometer length scale, can be dramatically different 
from the corresponding bulk material due to their very 
large surface-to-volume ratio. The unique optical and 
mechanical properties of such nanoclusters hold great 
promise for the development of nanotechnology. Gold 
nanoclusters, for example, have found application in such 
diverse fields as nano- lithography Q , catalysis Q , nano- 
bioelctronic devices |3|, and ion detection [4|. It is there- 
fore of great interest to study the behavior and stability 
of such nano-structures. 

In this work we consider gold nanorods of low aspect 
ratio ~ 3. Using molecular dynamic (MD) simulations of 
a rod with a few thousand atoms, we find that the sta- 
bility of the nanorod upon heating is governed primarily 
by the energetics of it's surface. When we start with a 
rod which has large {110} surface facets, the roughening 
of these facets upon heating leads to a shape transfor- 
mation of the rod to a shorter and wider structure. The 
surface reconstructs to form higher stability {111} facets, 
while the fee interior completely reorients to align with 
the new facet planes. When we start with a rod which 
has predominantly {111} facets on its surface, the rod 
remains stable until melting. 

The initial rod configuration that we consider, shown 
in Fig.^ is one which has been proposed || to apply 
to recent experiments 0, 13 on the laser heating of 
gold nanorods. The interior of the rod is a pure fee lat- 
tice. The surface of the rod consists of four large {100} 
and four large {110} facets oriented parallel to the rod 
axis. The ends of the rod have a {001} facet and four 
small {111} facets connecting the {110} and the {001} 
facets. These experiments found that, upon heating, such 
rods underwent a shape transformation to bent, twisted, 
shorter, wider, and ^-shaped clusters. Transmission elec- 
tron miscroscopy studies [8j observed point and planar in- 
ternal defects to accompany such shape transformations. 
Recent MD simulations Q of such rods, using a contin- 
uous heating procedure meant to model the laser heat- 



ing of experiments, found similar shape transformations. 
These simulations found the shape transformation to be 
accompanied by a structural change in which planes of 
interior atoms shift, converting local fee structure to hep. 
The extent and stability of these interior rearrangements 
was found to depend upon both the heating rate and the 
number of atoms in the cluster, but no specific mecha- 
nism or energetic argument for this structural rearrange- 
ment was proposed. In this paper we present new simula- 
tions carried out with a much slower "quasi-equilibrium" 
heating that allows the rod more time to approach its 
lowest free energy configuration. Our results make it 
clear that it is the energetics of the surface that is driving 
the shape and structural transformation. 

We use the empirical "glue" potential to model the 
many body interactions of the gold atoms in our simu- 
lated nanorod, and we integrate the classical equations 
of motion for the atoms using the velocity Verlet algo- 
rithm with a time step of 4.3 fs. However instead of 
increasing the kinetic energy at each MD step to model 
continuous heating as in Ref. , we now use the Gaussian 
isokinetic thermostat ^3] to keep the total kinetic energy 
fixed at a constant temperature T; after each MD step, 
all velocities are rescaled by a constant factor so as to 
keep ((l/2)mv 2 ) — (3/2)fesT fixed. Our procedure con- 
serves total linear and angular momentum, which are set 
to zero, so that our rod does not drift or rotate through- 
out our simulation. At each fixed T we carry out 10 7 MD 
steps, for a simulated time of 43 ns, before increasing the 
temperature in jumps of 100 K. Our effective heating 
rate is therefore ~ 2.3 x 10 9 K/s, more than three orders 
of magnitude slower than the continuous heating rates of 
2-7 x 10 12 K/s used in Ref. 0. We use a rod of N = 2624 
atoms with initial aspect ratio of 3, as shown in Fig.^ 
The length of the rod, parallel to its long axis, is 7.38 nm 
and its crossectional area has a diameter of 2.46 nm. We 
do a short equilibration for 430 ps (10 5 MD steps) at 5 K 
in order to relax the surface atoms of the rod from their 
initial fee positions, before starting to heat the rod. 
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FIG. 1: Initial configuration of the gold nanorod with N = 
2624 atoms and aspect ratio 3: (a) side view, (b) end view 
down the long axis, (c) cross-sectional view parallel to long 
axis. In (a) and (b), yellow atoms represent {111} facets, 
green atoms {100} facets, and red atoms {110} facets; white 
atoms are on the edges. In (c), yellow atoms are those with a 
local fee structure, white atoms are on the surface; the cross- 
sectional view in (c) shows atoms in the square arrangement 
of a {100} plane of the fee lattice. 

As a signature of the shape change of our nanoclus- 
ter we measure the radius of gyration, r g , defined by 
r g = (l/N) J2i \ r i — r c\ 2 , where is the position of atom 
i and r c is the center of mass. In Fig.|2t we plot our re- 
sults for r g as the system is heated; the blue curve is 
for our above "quasi-equilibrium" heating. The vertical 
dotted lines separate bins of constant temperature simu- 
lation, where the temperature is equal to the value at the 
left end point of the bin; the data plotted within each bin 
represents the instantaneous value of r g as a function of 
increasing time at the constant temperature. At the end 
of each bin the temperature is increased by a jump of 100 
K. We plot our data this way, instead of as an average 
value at each T, to highlight that significant shape relax- 
ation occurs even at constant T. For comparison, we plot 
r g for the continous heating (red curve) of Ref. |9j for the 
heating rate of 7 x 10 12 K/s. We see that the curves are 
qualitatively similar, with the onset of a plateau around 
800 K, however the present quasi-equilibrium heating al- 
lows the rod to relax to smaller r g values, before the rod 
melts at T ~ 1200 K. The decrease in the radius of gyra- 
tion reflects the shape transformation to a shorter wider 
rod of smaller aspect ratio. 

To investigate the local structure within the cluster, 
we use the method of bond orientational order parame- 
ters 01 . These parameters measure the orientation of 
bonds connecting a given atom to its nearest neighbors, 
and provide a convenient means of determining the lo- 
cal crystalline structure of an atom. In particular, we 
measure the 6-fold and 4- fold orientation parameters, Qe, 
Wq, Q4, and W4. We refer the reader to the literature for 
their definitions an d m Table[I]we give their values 

for several periodic three dimensional crystal structures. 
We also, in Table I, give the values of these parameters as 
computed for atoms on particular low index planar sur- 
faces of an fee bulk crystal; for these two dimensional pa- 
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FIG. 2: (a) Radius of gyration r g vs. temperature T for the 
rod of Fig.0 quasi-equlibrium heating (blue) compared to 
the continuous heating of Ref. 9] (red); continuous heating 
of the rod with structure of Fig.|3] (green) . (b) Bond order 
parameters, Q 4 , Qe, W4 and We, averaged over all atoms 
internal to the rod, vs. temperature T, for quasi-equlibrium 
heating of the rod of Fig.Q 



TABLE I: Bond order parameters for face-centered-cubic 
(fee), hexagonal close-packed (hep), simple cubic (sc), body- 
centered-cubic (bee), and low index fee planes. 



Geometry 


Q4 


Qe 


W 4 


W 6 


fee 


0.190 94 


0.574 52 


-0.159 32 


-0.01316 


hep 


0.097 22 


0.484 76 


0.13410 


-0.012 44 


sc 


0.763 76 


0.353 55 


0.159 32 


0.01316 


bee 


0.082 02 


0.500 83 


0.159 32 


0.01316 


{110} 


1 


1 


0.13410 


-0.093 06 


{100} 


0.829 16 


0.586 30 


0.124 97 


-0.007 21 


{111} 


0.375 00 


0.740 83 


0.13410 


-0.046 26 



rameters we average only over bonds connecting an atom 
with its neighbors in the specified plane. In Fig.l^b we 
plot these order parameters for our rod, averaging over 
only atoms internal to the rod (i.e. we exclude surface 
atoms since these have fewer nearest neighbor bonds). 
As for r gi we plot our data as the instantaneous value 
as a function of increasing simulation time, for bins of 
constant temperature (indicated by the dotted vertical 
lines). 

Comparing with Tabled we see that the rod maintains 
its fee structure until about 400 K. Then, from around 
400 K to about 800 K, there is a rise to positive values in 
W4, and a decrease in Q 6 and We, suggestive of a more 
hep-like structure. Above 800 K, the values return to 
their fee-like values. 

We now focus on the structure of the rod in the high 
temperature plateau region where r g stabilizes to a con- 
stant. In Fig. |3 we show the configuration of the rod at 
T = 900 K, in the constant plateau region before melt- 
ing. The views of the rod shown in Figs.[BK>b,c are the 
same orientations as shown for the initial configuration 
in Figs.^,b,c. In order to better illustrate the order 
of the rod, we first pick an instantaneous configuration 
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FIG. 4: Average cross-sectional shape, viewed down the long 
axis of the rod, for different temperatures. 



FIG. 3: Configuration of the nanorod after quasi-equilibrium 
heating to 900 K. (a) side view, and (b) end view down the 
long axis, after peeling away the surface and the first sub- 
surface layer; yellow atoms are {If 1} facets, green atoms are 
{f00} facets, and white are edge atoms, (c) cross-sectional 
view parallel to the long axis; yellow atoms have a local fee 
structure, green atoms have a local hep structure, and white 
atoms are neither. The cross-sectional view in (c) shows 
atoms in the close-packed hexagonal arrangement of a {111} 
plane of the fee lattice. 

sampled from the middle of the T — 900 K simulation, 
and use the conjugate gradient method [l^j to quench 
local thermal fluctuations. At such high temperatures, 
the surface can be partially disordered compared to the 
interior, due to the diffusion of atoms on and near facet 
edges and vertices We therefore use the cone al- 

gorithm |l6| to identify and peel away atoms on the sur- 
face and in the first sub layer below it, and in Figs.|3p,b 
show the configuration of the second sub layer of the rod. 
We see a very regular shape covered almost completely 
with stable {111} facets. Based on the values in Tabled 
we use the following criteria to identify atoms in this 
layer as belonging to particular low index planes: {111} 
if 0.7 < Q 6 < 0.9 and -0.08 < W 6 < -0.02; {100} if 
Qe < 0.7 and W e > -0.02; and {110} if Q 6 > 0.9 and 
Wq < —0.08. Atoms in Figs.EJi,b have been colored ac- 
cordingly. 

The cross-sectional view in Fig. [2: shows an almost 
pure fee interior, as was the case for the initial config- 
uration, however we now see a close packed hexagonal 
structure characteristic of a {111} plane of the fee lattice, 
rather than the {100} plane seen in the cross-sectional 
view of Fig.nj;. We thus see one of our main results: 
in order to align with the new {111} surface facets, the 
bulk fee structure has completely reconstructed itself to 
a new orientation. For interior atoms, we use the fol- 
lowing criteria to identify the local crystal structure: fee 
if Qa > 0.17 and W 4 < -0.10; hep if Q 4 < 0.13 and 
W4 > 0.07. Atoms in Figs.|St have been colored accord- 
ingly. 

To see how the rod evolves from its initial configuration 
(Fig.[TJ to its reconstructed shape (Fig.JSJ, we consider 



the average cross-sectional shape in a plane transverse to 
the long axis of the rod. We compute this average shape 
as follows. For each instantaneous configuration we first 
eliminate all atoms on the end caps of the rods, and all 
interior atoms of the rod, and then project the remain- 
ing surface atoms into the xy plane, perpendicular to the 
long axis of the rod. Placing the origin at the resulting 
center of mass, we divide the plane into 100 equal polar 
angles, and then compute the average position of all sur- 
face atoms in each angular division. This result is then 
averaged over 1000 different instantaneous configurations 
sampled uniformly throughout the simulated time of 43 
ns at each temperature T. We plot the resulting average 
cross-sectional shapes, for several different T, in Fig. 01 

At low T we see the octagonal cross-section of the ini- 
tially constructed rod of Fig.^ with the flat edges rep- 
resenting the initial {100} and {110} facets. The shape 
stays roughly the same until about 400 K. Somewhere be- 
tween 300 — 400 K, the shape becomes rounder and the 
initial flat edges disappear. As T increases further, the 
cross-sectional area grows, representing the shape trans- 
formation to a shorter and wider rod of lower aspect ratio, 
and we see new flat facets develop and grow in new di- 
rections. At 900 K we see the fully faceted shape shown 
in Fig.E2 

The disappearance of the initial flat facets, with the 
resulting rounding of the average shape, is a signature 
of the roughening transition of those surfaces [TjJ. For 
macroscopic gold samples, it is known from experiments 
that the {110} surface roughens at 680 K 0, while the 
{100} surface disorders at ~ 1170 K [ljj, below the bulk 
melting temperature of 1337 K. In contrast, the {111} 
surface is believed to remain stable up to, and even above, 
the bulk melting ji^J] ■ The initiation of the shape change 
that we find in our nanorod at 400 K is most likely a con- 
sequence of the roughening transition of the {110} facets, 
which has been shifted to lower temperature due to large 
finite size effects in our relatively small rod (just as the 
bulk meltin g tr ansition can be greatly reduced by finite 
size effects 2 1] ) . After this roughening, the surface re- 
constructs to form mostly lower free energy {111} facets, 
which remain stable until melting. Note that the fully 
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FIG. 5: Cross-sectional view for various temperatures. Atoms 
are colored according to local crystal structure: fee is yellow, 
hep is green, neither is white. 



formation, we study the stability of a gold nanorod with 
an aspect ratio of ~ 3, but with an initial structure simi- 
lar to that of Fig.|21 with a surface predominantly covered 
by {111} facets. We use a continuous heating MD simu- 
lation with a heating rate of 7 x 10 12 K/s to model laser 
heating experiments, for a rod with 3411 atoms. Our 
results for r g vs T are shown in Fig.|2Jt (green curve). 
Unlike the intial rod of Fig.EJ we now find that the rod 
remains stable, with no significant shape or structural re- 
arrangement, up until the rod melting temperature. We 
conclude that the stability of gold, and presumably other 
metallic, nanorods is crucially dependent upon the struc- 
ture of the rod surface. 

This work was funded in part by DOE grant DE-FG02- 
89ER14017. 



faceted cross-sectional shape at 900 K contains four large 
sides and two short sides; the former are the {111} facets, 
while the latter are {100} facets. By 1100 K, these {100} 
facets have been replaced by a smoothly curved surface. 
We infer that this is due to the disordering transition of 
the {100} surface, reduced somewhat in temperature due 
to finite size effects. 

In order to see how the interior fee structure of the 
rod reconstructs itself to a new orientation, we show in 
Fig.|S]cross-sectional views of the rod at various tempera- 
tures. We color the atoms according to their local crystal 
structure, using the criteria given above: fee is yellow, 
hep is green, neither is white. Initially, the interior is 
pure fee, oriented so that the cross-sectional view shows 
a {100} plane of atoms. As temperature increases, we see 
that the shape and structural transformation is accom- 
panied by the appearance of hep planes inside the rod 
interior, due to the sliding of {111} planes. As tempera- 
ture further increases, the surface becomes less ordered, 
and more {111} planes with different orientation slide. 
Around 800 K, the surface has reordered and the inte- 
rior fee lattice has reoriented so that the cross-sectional 
view now shows a predominantly {111} plane of atoms. 
At 1100 K, the interior has completely reordered to pure 
fee, but with the new orientation. 

Such behavior as described above may well exist in 
other simple elemental metals. We note that many such 
metals similarly have a roughening transition Tr for the 
{110} surface that is significantly below the bulk melting 
T m . Silver, for example, has Tr = 600 K and T m = 1235 
K, with a similar ratio of T^/T m as gold 22]. Lead has 
Tr = 415 K and Tm = 601 K, for a somewhat larger 
T R /T m than gold 0. 

Finally, in order to verify that the roughening of the 
{110} facets, rather than just the minimization of total 
surface area, is indeed the mechanism for the shape trans- 
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